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The selection of basic variables in current-density functional theory and formal properties of the 
resulting formulations are critically examined. Focus is placed on the extent to which the Hohenberg- 
Kohn theorem, constrained-search approach and Lieb's formulation (in terms of convex and concave 
conjugation) of standard density-functional theory can be generalized to provide foundations for 
current-density functional theory. For the well-known case with the gauge-dependent paramagnetic 
current density as a basic variable, we find that the resulting total energy functional is not concave. 
It is shown that a simple redefinition of the scalar potential restores concavity and enables the 
application of convex analysis and convex/concave conjugation. As a result, the solution sets arising 
in potential-optimization problems can be given a simple characterization. We also review attempts 

Oto establish theories with the physical current density as a basic variable. Despite the appealing 
physical motivation behind this choice of basic variables, we find that the mathematical foundations 
QQ °f the theories proposed to date are unsatisfactory. Moreover, the analogy to standard density- 

functional theory is substantially weaker as neither the constrained-search approach nor the convex 
i i analysis framework carry over to a theory making use of the physical current density. 
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I. INTRODUCTION 



Density-functional theory (DFT) constitutes one of the most popular methods in quantum chemistry. The foun- 
dations of DFT rest in particular on three contributions: First, the Hohenbcrg-Kohn (HK) theorems established 
a one-to-one mapping between a set of scalar potentials and a set of ground-state densities as well as a variation 
principle based on the density [1] . Here, the density is the charge density (strictly, the negative of the charge density 
in units of the elementary electron charge). Second, the Levy-Lieb constrained-search expression provided a formal 
but explicit expression for the intrinsic energy (the universal density functional) and clarified significant fundamental 
points [2 . Third, Lieb further generalized the universal functional to a convex functional represented in terms of a 
Legendre-Fenchel transform 3J. From a mathematical point of view, Lieb's formulation is particularly attractive as 
it allows application of convex analysis to establish several properties of the intrinsic energy functional [4-6 . Ad- 
ditionally, Lieb's framework has made feasible practical calculations of approximations to the exact intrinsic energy 
functional and adiabatic connection curves [71-11 1| , enabling detailed comparisons of the properties of approximate and 
near-exact density functionals to be made [T2HT7] . 

Standard DFT, involving universal energy functionals of only the charge density, is limited to the treatment of 
physical systems that may be represented as eigenstates of Hamiltonians that differ only in their scalar potentials. 
To treat systems subject to an external magnetic field, it is necessary to introduce an additional dependence on the 
magnetic field or its associated vector potential into the Hamiltonian. Consequently, a dependence on a corresponding 
variable apart from charge density is needed in the universal energy functional. In magnetic-field density functional 
theory (B-DFT) this is resolved by constructing a family of density functionals — one for each external magnetic 

field pa eng. 

In the present work, we consider the alternative current density- functional theory (CDFT), where the additional 
variable is either the paramagnetic current density or the physical current density. We restrict our attention to 
non-relativistic formulations and most of the discussion will for simplicity not be concerned with densities or density- 
contributions arising from spin-degrees of freedom. We term the variables on which the energy functionals explicitly 
depend the basic variables and make a distinction between basic densities and basic potentials. Many choices of basic 
densities are conceivable [201 HI]; we- require only that the choices result in useful density- functional theories. Our 
perspective thus differs from that in recent works on CDFT by Pan and Sahni [2"2T|25| . who restrict the term basic 
variable to variables that admit an HK theorem. Although it appears naturally in the generic framework introduced 
by Ayers and Fuentealba [21], the possibility of choosing basic potentials other than the standard electromagnetic 
potentials and fields has not previously been explored in detail. 

By far the most developed form of CDFT is that due to Vignale and Rasolt [26l [27] . who use the charge and 
paramagnetic current densities as basic variables. For these variables, a Kohn-Sham approach has been formulated |26j 
with an associated adiabatic-connection [55], virial and scaling relations [2"8H5T] analogous to standard Kohn-Sham 
DFT. In addition, optimized-effective-potential (OEP) approaches based on this formulation of CDFT have been 
presented to treat non-collinear magnetism 32 - 35j and extensions to time-dependent CDFT have been considered [3BT - 

m\. 

However, in CDFT based on the charge and paramagnetic current densities as basic variables, no HK-type theorem 
exists and the consequences of this have been extensively discussed in the literature [43]. In the present work, we 
examine this question for CDFT in some detail, demonstrating how convex analysis of the underlying universal density 
functional can be a significant aid in clarifying the relationship between basic variables of CDFT and the potentials. 
A CDFT featuring the gauge-invariant physical current density (rather than the paramagnetic current density) as a 
basic variable is appealing from a physical perspective and is therefore also considered here. Specifically, we examine 
the formulations due to Diener jH] and Pan and Sahni [2"!3] . 

We begin in Sec.|ll]by introducing notation related to sets of basic potentials, basic densities, and mappings between 
them. In Sec. [Hl[ we consider CDFTs that use the charge and paramagnetic current densities as basic variables — in 



particular, Sec. HI E| establishes the concavity of a universal density functional based on these variables and Sec. Ill H 



outlines the opportunities that this formulation affords for numerical studies of this functional. Next, in Sec. |IV[ the 
use of the charge and physical current densities as basic variables is considered and two previous formulations [22|, [44] 
are examined. Our concluding remarks are presented in Sec. [V] 



II. A REVIEW OF DFT 



Before discussing CDFT, we briefly review standard DFT, with emphasis on Lieb's treatment based on convex 
conjugation [3J. The concepts and techniques of convex analysis introduced here are well suited to the study of DFT 
and will later be used in our discussion of CDFT. Some background is also given in the Appendix. 
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We consider a system of N electrons with an electronic Hamiltonian of the form (in atomic units) 

k k 

where p/j = — iVfc is the canonical momentum operator of electron fc, v(r) is the external potential at position r, and 
W = ^2 k<t is the two-electron Coulomb repulsion operator. The state of the system is described by a density 
matrix T, which is a convex combination of normalized iV-electron pure-state density matrices 

r = Y / MA)(A\, \>o, ]Ta 1 = i, (2) 

i i 

where the wave functions tpi are antisymmetric in the space and spin coordinates = (r^Ok) of the N electrons. 
The electron density associated with such density matrices is given by 

p(r) = X>ft(r), Mri) = N U*^dr^, (3) 

where the volume element is dr_i = dcridx2 • • -dxjv, i.e., the integration is over all N spin and spatial coordinates 
except ri. The ground-state energy is obtained from the Rayleigh-Ritz variation principle, 

E[v}= MTt(TH[v}), (4) 

where the minimization is over all ./V-electron density matrices. 

An infimum rather than a minimum is taken in Eq. Q since v may or may not support an ./V-electron ground 
state. The set of potentials that support one or more iV-electron ground states (and for which therefore the infimum 
is attained) is denoted by Vat; the potentials in Vn are sometimes said to be p-representable. Conversely, a density 
that is an ensemble ground-state density for some potential v € Vat is said to be (ensemble) ^-representable; the set 
of u-representable densities is denoted by Bn- For convenience, we shall also refer to p-representable potentials and 
u-representable densities as ground-state potentials and densities, respectively. 

In the constrained-search formalism of DFT, we write the Rayleigh-Ritz variation principle as an HK variation 
principle, 

E[v] = inf (F[p] + (p\v)) , (5) 

where In is the set of iV-representable densities — that is, the set of the nonnegative densities p with J p(r)dr = N 
and with a finite von Weizsacker kinetic energy. The Lieb constrained-search functional F is given by 

F[p] = inf Tr(r#[0]), (6) 
rn-p 

where the notation r i— > p indicates that the minimization is restricted to density matrices T that reproduce the 
density p. If p is not 7V-representable, no such T exists, and F[p) = +oo by definition. The HK variation principle 
in Eq. ^ is well defined for all potentials v that have a finite pairing with every p £ Xjv, 

(p\v) = Jp(r) u(r)dr. (7) 

An especially attractive formulation of DFT is Lieb's formulation in terms of Legendre-Fenchel transformations or 
convex conjugation. This formulation is not only elegant, but also fits naturally in the well-developed mathematical 
field of convex analysis, allowing application of deep results of convex analysis to DFT, of which we will give some 
examples. 

Lieb's formulation of DFT begins with the observation that the ground-state energy E[v] is upper semi-continuous 
and concave in v and therefore may be represented by its conjugate function: Lieb's universal density functional F[p\. 
The ground-state energy and density functionals are then related as 

E[v] - inf [F[p] + (p\v)} , (8a) 
ps x 

F[p] = sup [E[v] - (p\v)} . (8b) 

vex* 
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Here, X is a Banach space (a complete normed vector space) that contains Ijy, and X* is its dual — that is, the set 
of all bounded linear functionals on X (also a Banach space), thereby ensuring that |(p|w)| < +00. Lieb identified 
X = L 1 n L 3 C In an d X* — L?l 2 + L°°, which contains, among others, all Coulomb potentials. The observations 
that E[v] and F[p] are upper and lower semicontinuous concave and convex functions, respectively, together with 
the identification of the Banach spaces X and X* are the key elements that place DFT within the setting of convex 
analysis. 

The duality of E and F apparent in Eq. ([8]) means that the same information is contained in either functional but 
encoded in different ways; this duality is emphasized by referring to E and F as the extrinsic and intrinsic energies, 
respectively, of the electronic system. The infimum and supremum expressions on the right-hand sides of Eqs 



and (|8bj) feature linear pairings of densities and potentials and are therefore by construction concave and convex, 
respectively. Indeed, a necessary and sufficient condition for E[v] (F[p\) to be upper (lower) semicontinuous and 
concave (convex) is the existence of such expressions [IS] . 

Unlike the Lieb density-matrix constrained-search functional, the Levy-Lieb constrained-search functional, defined 
in terms of pure states rather than density matrices, is not convex and hence not identical to the Lieb functional. 

In the present paper, we generalize Lieb's formulation of DFT to CDFT. In particular, we discuss the Vignale- 
Rasolt constrained-search functional as a generalization of the Lieb functional to systems in the presence of a vector 
potential. We shall see that such a generalization is possible after a redefinition of the scalar potential. However, we 
here leave aside technical questions such as lower semicontinuity and infimum and supremum domains (i.e., Banach 
spaces), which is the subject of future work; for a discussion of such mathematical issues within standard DFT, 
see Lieb [3] and Eschrig By contrast, the concavity of E and the convexity of F are essential properties for 
our discussion of CDFT. If E[v] happens to be non-concave, it cannot be represented by an expression like that in 
Eq. (8a I, even by allowing F to be non-convex. Therefore, no universal functional with a linear potential pairing can 
exist for non-concave energies (such as those of excited states of the same symmetry as the ground state). 



III. THE PARAMAGNETIC CURRENT DENSITY AS A BASIC VARIABLE 



A CDFT with the paramagnetic current as a basic density was considered in the seminal work of Vignale and 
Rasolt [26l [27] . In their formulation of CDFT, the basic potentials are the standard electromagnetic potentials (v, A) 
and the basic densities are the charge density and paramagnetic current density (p,j P ). We shall here first review 
their theory and then discuss an alternative formalism, based on a redefinition of the basic scalar potential. 



A. Preliminaries 



We consider electrons subject to time-independent external electromagnetic fields E(r) = — Vw(r) and B(r) = 
V x A(r), represented by the scalar potential v(r) and the vector potential A(r), respectively. For potentials (u, A), 
we introduce the equivalence relation 

(v', A') ~ (v, A) (VV, V x A') = (Vt>, V x A), (9) 

which defines equivalence classes of potentials that differ only by a static gauge transformation, thereby representing 
the same external fields. 

We note that a general gauge transformation of v and A is given by v' — v — df jdt and A' = A + V/, for some 
arbitrary gauge function f(r,t). If A is to remain static after the transformation, we must require that f(r, t) — 
x( r ) ~ ct, where c is constant. It follows that a general time-independent gauge transformation is given by v' = v + c 
and A' = A + V%, where the constant c and the function x(r) are independent. Therefore, the equivalence relation 
in Eq. ([9| holds if and only if there exists a constant c and a sufficiently well-behaved gauge function %(r) such that 
v' = v + c and A' = A + V\- 

In the presence of a vector potential, the electronic Hamiltonian in Eq. ([I]) is modified by replacing the canonical 
momentum operator p& = — iV^ by the mechanical (kinetic) momentum operator = — iVfc -I- A(rfc), yielding 

H[v,A] = ±J2rt+T,< r >*) + W - ( 10 ) 

k k 

We have here omitted the spin-dependent term, X)fcB( r <0 ■ S, from the Hamiltonian. By analogy with Eq. Q, the 
Rayleigh-Ritz variation principle in the presence of a vector potential is given by 



E[v,A] = inf Tr(ri/[u, A]), 



(11) 
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where the minimization is over density matrices T containing N electrons, see Eq. ([2]). An infimum rather than a 
minimum is taken to ensure that the energy is well defined also when (v, A) does not support a ground state. We 
denote the set of all potentials (v, A) that support a ground state with this Hamiltonian by 

V N = {(v,A)\H[v,A] has a g.s.} (12) 



and also introduce the related set 



U N = {{v+ \A 2 ,A)\H[v,A] has a g.s.} (13) 

in preparation of a reparameterization of the scalar potential that will be introduced later. 

In the presence of a vector potential, the ensemble ground-state charge densities p are as before given by Eq. (J3j) . 
Regarding the induced currents, we distinguish between the paramagnetic current density and the physical current 
density. The former is defined as 

j p (n) = Re^Ai J ip*p k iPi dr_ x . (14) 

ik 

The paramagnetic current density is gauge-dependent and unobservable. The physical current density is given by 

j(ri) =Rey^Aj / tp*7r k ipidT-i, (15) 

ik 

and satisfies the relation j = j p + pA. Unlike the paramagnetic current, the physical current is gauge-invariant. 
Finally, the sets of paramagnetic and physical u-representable ground-state densities are denoted by 

B P n = {(pJp)I(pJp) is S- s - den - of some H[v, A]}, (16) 
#;v = {(P,j)|(p,j) is g-s. den. of some H[v,A]}, (17) 

where both mixed and pure states are allowed. 



B. Do paramagnetic densities determine potentials? 



The HK theorem of standard DFT states that the ground-state density p determines the scalar potential v up 
to a constant shift. Hence, two potentials that differ by more than a constant shift cannot give rise to the same 
ground-state density. This fact establishes a mapping from ground-state densities to potentials. 

Vignale and Rasolt established that two different potentials (vi,Ai) ^ (u 2 ,A 2 ) with different ground-state wave 
functions ipi ^ ip2 cannot give rise to the same paramagnetic ground-state density (/9,j p ). However, this does not 
establish an analogue of the HK theorem for CDFT [43] since different potentials Ai) ^ (u 2 , A 2 ) can map to the 
same density in B P N via the same wave function, *f?[vx, A x ] = *f?[v2, A 2 ]. Vignale and Rasolt's result can be extended 
to a form that applies also when the ground states are degenerate: Let ipi be a ground state of H\v-y, AJ and let ip2 
be a ground state of H[vz, A 2 ] . If ip\ and tp2 give rise to the same paramagnetic density, (/tyi , jpj-^i ) = (p^a; jp;^)* 
then ip2 is also a ground state of H[vi, A\] and ipi is also a ground state o/iJ[w 2 , A 2 ]. 

The above statement appears to be as close as one can get to a HK-like result for paramagnetic densities (/0,j p ). 
A CDFT formulated in terms of the paramagnetic current density thus cannot be based on a formal mapping from 
ground-state densities to potentials. On the other hand, rigorous formulations of standard DFT do not rely on the 
HK mapping from densities to potentials. These rigorous formulations can be extended to the paramagnetic current 
density as a basic variable. The absence of an HK-type theorem in CDFT is therefore not a serious impediment. 



C. The standard electromagnetic potentials as basic potentials 



An important property of E[v] in Eq. (8a) is its concavity in v, which established duality with the universal density 
functional F[p\. However, unlike E[v], the energy functional E[v, A] in Eq. (11) is not concave. The non-concavity 
of E[v,A] is apparent in, for example, any diamagnetic ground state at vanishing external magnetic field. Such a 
ground state has a negative definite magnetizability tensor x and, when restricted to weak uniform magnetic fields 
B sa 0, the energy is a convex function Eq — ^B T xB in B and therefore in A. 



G 



In more detail, consider a one-electron system confined to the (two-dimensional) xy-plane, subject both to a uniform 
magnetic field along the z-axis and to a harmonic-oscillator potential. Parameterizing the scalar and vector potentials 
under consideration as 

v HO (r;k) = l -k{x 2 +y 2 ) 1 A ± (r; B) = ^Be z x r, (18) 

we obtain the following Hamiltonian 

H[vh , AJ = \ P 2 + X -BL Z + i7 HO (r; k) + \a±(t; B) 2 

= \p 2 + \ BLz + I (4k + B 2 ) (x 2 +y 2 ), 

where L z is a good quantum number. For k > k m - m > and some finite interval \B\ < B max , the ground state has 
L z = and application of a magnetic field has exactly the same effect as the introduction of a harmonic-oscillator 
potential. For these potentials, the ground-state energy is 

E[v HO , Ax] = E[v m + \A\, 0] = y/k + ±B*. (20) 

Note that the right-hand side is concave in k > and convex in B. Hence, on the restricted set of potentials spanned 
by k > fc m ; n and \B\ < B maK , the functional E[vao, Aj_] is not only non-concave but convex in its second argument. 
On a larger domain, the functional is neither concave nor convex. We conclude that E[v, A] cannot be represented 
by a conjugate functional in the manner of Eq. ([8]). However, this does not preclude a constrained-search formulation 
of CDFT, as discussed in the next subsection. 



D. CDFT by constrained search 



Rewriting the Rayleigh-Ritz variation principle in Eq. (11) by analogy with the constrained-search approach of 
standard DFT in Eq. we obtain an HK-type variation principle for a system in the presence of a scalar and vector 
potential, 

E[v, A] = inf [FvrIp, j p ] + (p\v + \A 2 ) + (j p |A)] , (21) 

P dp 

where the Vignale-Rasolt constrained- search functional is given by 

FvK[p,j P ]= inf Tr(r(ip 2 + W0) (22) 
ri -HPJp) 

and we have introduced the following notation for the pairing between a current density and a vector potential: 

(j P |A) = J j p (r)-A(r)dr. (23) 



Like the Lieb constrained-search functional given in Eq. (22), the Vignale-Rasolt constrained-search functional 
in Eq. (22) is universal in the sense that it does not depend on the potential (v,A), only on the density (/0,j p ). 



Another important characterization of the Lieb functional is its convexity. To examine the convexity of Fys,, let 
G°ijjpi) an d (P2)jp2) be arbitrary (Lebesque integrable) functions and let < A < 1, \i = 1 — A. We then obtain 

Fvr[^Pi + MP2, Ajpi + ^j p2 ] 

< inf. Tr((Ar 1+Ai r 2 )(ip 2 + W0) 

rin-(pi j P i) (,^4J 
r 2 >-Kp2 j P 2) 

= AF V r[pi, jpl] + /J,Fy R [p2, j p 2], 

demonstrating that J^vr is convex in (p,j p ). The key point in establishing the inequality above is to restrict the 
infimum over all density matrices T h-> (p,j p ) to an infimum over all matrices of the form T = XTi + /ir 2 where 
Ti i-> (pi,jpi) and T 2 i-> (/0 2 ,j p2 ), thereby overestimating the infimum. 

Given that the Vignale-Rasolt functional is convex, it is uniquely represented by a convex conjugate functional. 
For the Lieb functional F[p], the conjugate is the concave ground-state energy E[v}. However, since E[v, A] is not 
concave, it cannot be the conjugate to -FvR[/9,j p ]. In the following, we identify the energy conjugate to Fvr[p,J p ], 
thereby arriving at a Legendre-Fenchel formulation of CDFT. 
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E. CDFT by convex conjugation 



Inspection of the general expression in Eq. (21 1 and the harmonic-oscillator example suggest the introduction of a 
new basic scalar potential, 

u = v + \A 2 . (25) 

The choice of (it, A) as basic potentials and (p,j P ) as basic densities results in a theory where the HK variation 
principle takes the form of a Legendre-Fenchel transformation with a linear pairing (p\u) + (j p |A) of the densities 
and potentials: 

E[u, A] = inf Tt(TH[u, A]) 

= inf [F VR [p, j p ] + (p\u) + (j p |A)] . (26) 



Here we have introduced the notation 



H[u,A] = H[u- \A\ A], (27) 



2 

E[u,A] = E[u- \A 2 ,A]. (28) 

The energy E[u, A] is now by construction concave, allowing it to be generated from the convex intrinsic energy 
Fvr[p, j p ] by a reverse Legendre-Fenchel transformation: 

Fvr[ P , j P ] = sup [E[u, A] - (p\u) - (j p |A)] . (29) 

u.A 

Thus, by a change of variables from v to u = v + \A 2 ; , we have restored the conjugate relation between the extrinsic 
and intrinsic energies characteristic of standard DFT. 

Strictly speaking, for i*vR.[/>, j p ] and E[u, A] to form a conjugate pair, we must specify their domains in the form of a 
Banach space X p and its dual X*, respectively. Furthermore, we must demonstrate lower and upper semi-continuity 
of the intrinsic and extrinsic energies, respectively. However, regarding the domains, we note here that, in Lieb's 
formulation of standard DFT, the vector space X* = L 3 / 2 + L°° of potentials does not contain every potential with 
a square integrable ground state (e.g., it does not contain harmonic potentials), but X* does contain all Coulomb 
potentials, thus covering most systems of interest. A similar compromise is expected for CDFT: we cannot expect to 
identify a vector space X p of densities such that its dual X* includes all ground-state potentials Vjv- In particular, 
since gauge transformations may produce potentials A + Vx that are arbitrarily ill-behaved at infinity, we expect 
some gauge restriction to be necessary. 



F. Subdifferentiability in CDFT 

In general, the functionals E and Fvr are not differentiable. Consequently, we cannot characterize the ground- 



state densities in the CDFT HK variation principle of Eq. (26) in terms of functional derivatives. On the other 
hand, in convex analysis, the proper tool for characterizing minimizers and maximizers are sub- and supergradients, 
respectively. We here introduce sub- and supergradients in the context of CDFT. 



We begin by noting that an immediate consequence of the CDFT variation principles in Eqs. (26) and (29) is 
Fenchel's inequality, 

E[u, A] < Fynip, j p ] + (p\u) + (j p |A), (30) 

valid for any choice of potential (v, A) € X* and density (p,j p ) € X p . Moreover, equality holds if and only if 
(/°Jp) € &n 1S a ground-state density belonging to (u, A). To characterize ground-state densities and their potentials 
mathematically, we use the concepts of subgradients and subdifferentials. A subgradient of FyR at (po>jpo) S AT p is 
an external potential — (uq, A ) € X* for which the inequality 

^Vr[pJ p ] > ^vr[PoJ p o] - (p-PoK) ~ (j P -j P o|A ) (31) 
holds for all (p,j p ) € X p ; see Fig.JT] Clearly, all potentials — (u , A ) for which the density (poJpo) is a minimizer in 



Eq. (26) are subgradients of fyR at (pojjpo)- The set of all subgradients at (pcbjpo) is known as the subdifferential of 




FIG. 1. Schematic illustration of the subdifferential 8Fvn[po, J p o]' The set of all slopes (potentials — (it, A)) of just-touching 
tangent planes entirely below the graph at (po,jpo). One particular tangent, with slope — (ito, Ao), is shown, while others are 
indicated by dashed lines. A similar illustration can be made of the superdifferential BE[u, A]. 



Fvr at (poJpo) and is denoted by 9-FVr[po Jpo] C X*. Hence, to within a minus sign, the subdifferential at (poJpo) 
is the collection of all external potentials that have the same ground-state density (p , j p o). 

Analogously, we consider the concave energy functional E[u, A] and its supergradients. In general, (po,j p o) £ X p 
is a supergradient of E at (it , A ) C X* if and only if — (poJpo) is a subgradient of the convex functional —E at 
(i*o,Ao). Hence, (po,jpo) is a supergradient of E at (uo, Ao) if the inequality 

E[u, A] < E[u Q , A ] + (po\u - u ) + (j p0 |A - A ) (32) 

holds for all (u,A) G X*. This condition is satisfied precisely when (poJpo) is the density arising from a (possibly 
degenerate) ground state of H[uq, Ao]. The superdifferential 8E[uq, Ao] is the collection of all supergradients of E at 
(u 0) A ). 

For all ground-state densities (pojjpo) and associated potentials (uq,Aq), we now have the following stationary 



conditions of the HK and Lieb variation principles in Eqs. (26 1 and (29), respectively: 



-(u , A ) € dFynlpo, j p0 ], (p , j P o) eB p N nX, (33) 
(po,jpo) e025[uo,Ao], {u ,A ) eU N DX*. (34) 

Importantly, these conditions are equivalent: — (uq, Aq) is a subgradient of -FVr at (po, jpo) if and only if (po, j p o) is a 
supergradient of E at (uq, Aq). Hence, instead of a one-to-one mapping between individual potentials and individual 
ground-state densities, the convexity and concavity of the intrinsic and extrinsic energies, respectively, establish a 
mapping between the convex sets {(p,j P )} C X p of degenerate ground-state densities and convex sets {(u, A)} C X* 
of potentials that give rise to identical ground-state densities; see Fig. [2] 

The sub- and superdiffcrentials are empty when no minimizer and maximizer exist in the corresponding optimization 



problems in Eqs. ( 26 ) and ( 29 1 . However, it is a general result of convex analysis that the subgradients (supergradients) 
of a convex (concave) function (under certain semicontinuity conditions) exist at a dense subset of the domain of the 
function. In CDFT, this result implies that the set of ground-state densities (po? jp.o) is dense in the set of all densities 
and that the set of potentials (mq, Aq) that support a ground state is dense in the set of potentials. 



G. Degeneracies in CDFT 



Consider now the case of degenerate ground-state densities in Eq. ( [26] ) . For a finite degeneracy Gd, the superdif- 
ferential of the ground-state energy in the external potential (i*o, Ao) is then a (Gd — l)-dimensional simplex with Gd 
pure-state densities (p , j p0 ) at the vertices: 

BB[ Uo , Aq] = co{( Po , j P ,o), |* = 1, G d }, (35) 
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FIG. 2. Illustration of considerations in Sec. |III F| Instead of a one-to-one mapping of individual potential pairs and corre- 
sponding ground-state density pairs, the sub- and superdifferentials of Fvr and E, respectively, maps (ground-state) densities 
into convex sets of (ground-state) potentials and potentials into convex sets of densities. 




FIG. 3. Illustration of the set dE[u, A] of ground state densities for potentials (u, A) with degeneracy Gd = 4. This is a 
convex set, a three-dimensional simplex, with vertices at (poJp,o)> i- e -i dE[u, A] = co{(pg, j P )}. In this example, we embed 
the simplex in R 3 , but in reality it is a three-dimensional subset of the density space X p . The corresponding set of potentials 
is also a convex set, but usually with a more complicated structure. See also Fig. [2] 



see Fig.p^for an illustration. Consequently, each ground-state density (po,j p o) may be written as convex combination 
of the Gd pure-state densities, 

Gd Gd 

(po Jpo) = E X *(Po 4o), = 1, A< > 0. (36) 

i=l i=l 

As is well known, such degeneracies are either accidental or caused by symmetries of the Hamiltonian. 

Consider next a ground-state density (po,jpo) with several maximizing potentials in the Lieb variation principle 
in Eq. (29 1. Like the superdifferential of the ground-state energy, the subdifferential of the Vignale-Rasolt density 
functional is a convex set. Let now (uq, Ao) be a potential with ground-state density (pojjpo) an d consider the family 
of external potentials 

MA), A (A)) = K>, A ) + X(Au , AA ), (37) 
for some potential (Au , AA ) and A € K. The associated Hamiltonians are given by 

H[u (X),A Q (X)}^H[u ,A }+XK, (38) 
K = ^[Auote) + \{\>i, AAofc)}]. (39) 

where {A, fl} denotes the anti-commutator. Clearly, if H[uo, Ao] and K commute, 

[B[u ,A ],K] =0, (40) 
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then H[uq(X), Aq(A)] has the same eigenstates for all A. In particular, as A changes continuously from 0, the external 



potentials in Eq. (37 1 will have the same ground-state density (poJpo) until a level crossing with the ground state 



occurs for H[u (XJ^A (X)}, see Capelle and Vignale [35] . We note that the commutation condition in Eq. (401 is 
sufficient but not necessary for the existence of degenerate maximizing potentials. 

There may be several, independent perturbing potentials that commute with the reference Hamiltonian H [uq, A ]. 
Assuming a finite degeneracy G p of these potentials, we may write the subdifferential of the Vignale-Rasolt density 
functional as a (G p — l)-dimensional simplex 

&Fvr[p 0i Jpo] = - co{(4, A*), \i = l, G p }. (41) 
All potentials with the ground-state density (po Jpo) ma Y then be written as a convex combination, 

(uo, A ) = Xi(ul A p0 ), £ Xi = 1, Xi > 0, (42) 

of the G p vertex potentials (u l , A l ). 

Consider now the special case of DFT, for which the HK and Lieb variation principles are given by Eqs. 



and (8b), where E[v] and F[p] depend only on the scalar potential and the electron density, respectively. Let vq be a 
potential with ground-state density pq. Clearly, the only scalar potentials v(r) that commute with the Hamiltonian 
H[vq] are the constant potentials c. It follows that the subdifferential of F at po is the convex set 

dF[p Q ] =-{«o(r) + c|ceM}, (43) 

which may be regarded as a one-dimensional simplex with vertices — vq(t) ± oo. In DFT, therefore, the ground-state 
density po determines the external potential «o( r ) uniquely up to an additive constant c, in accordance with the HK 
theorem. 

Returning to CDFT, consider next two potentials [u\, Ai) and (u 2 , A2) with the same ground-state density (po, j p o)- 
By the convexity of the subgradicnt of Fvr> all convex combinations 

(u, A) = A(«i, At) + (1 - A)(U2, A 2 ), < A < 1 (44) 

then have the same ground-state density. Recalling that U\ = Vt + \A\ , u 2 = v 2 + \A^, and u = v + \A 2 , the 
characterization of the non-uniqueness given in Eq. ( |44[ ) can be expressed in terms of the ordinary scalar potential v. 
If (vt, Ai) and (v 2 , A 2 ) give rise to the same density, then so do all potentials of the form 



v\ _ /Ai;i + (l-A)u 2 + |A(l-A)|Aa 
AI { AA X + (1-A)A 2 



(45) 



with < A < 1. However, this set is not a convex set and not a subdifferential, due to the use of the (v, A) rather 
than (it, A) variables. 

An advantage of the formulation of stationary conditions in CDFT in terms of sub- and superdifferentials is that 
differentiability is not required. In general, a sufficient condition for differentiability of a function at a point is that 
the function is continuous at this point and has a single sub- or supergradient there; in the absence of continuity, 
differentiability is not guaranteed. The ground-state energy is differentiable at all potentials (u, A) that have a 
nondegenerate ground-state density, whereas the Vignale-Rasolt functional is in principle nowhere differentiable since 
we may always add a constant term c to the potential without affecting the ground-state densities. However, assuming 
that this is the only cause of nondifferentiability of the potentials, we may in the absence of other degeneracies write 

Mr) ( ' ' £j P (r) AW ' l4bj 

and 

SE[u, A] SE[u,A] 

5u(r) ~ P( h SA(r) ~ Jpl h 10 



where (p,j P ) is the (non-degenerate) ground-state density of the potential (w, A). 
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H. Numerical calculation of Fvu[p,] p ] 



Besides allowing formal application of theorems in convex analysis to CDFT, the convex formulation given above 
has practical value. After linear programming and optimization of quadratic functions, optimization of convex and 
concave functions is the mathematically most well-characterized type of optimization problem. The fact that convex 
and concave optimization problems have a unique global optimum (either in the form a single point or a convex set 
of optimal points) and no additional local optima is of great value when devising practical optimization methods. 

In standard DFT, Lieb's formulation [3] of F[p) in terms of the Legendre-Fenchel transform of the concave energy 
functional E[v) has proven useful in the study of functionals of interest in Kohn-Sham theory [T2TU7| . In particular, 
the modulation of the two-electron interaction operator W by a parameter A such that, 

W\ = '^2,w\{r i j), w (rij) = 0, w^nj) = l/r#, (48) 

i<j 

allows us to represent the ground-state energy E\[v] at interaction strength A in terms of its conjugate functional 
F\[p\. The standard choice is w x (rij) = X/rij, but other choices are possible. If the density supplied to F\{p] is held 
fixed at a physical density generated by an appropriate (high-level, systematically improvable) ab-initio quantum- 
chemical methodology and if the value of the interaction-strength parameter is varied from to 1, then the adiabatic 
connection [TKTTj between the Kohn-Sham (A = 0) and physical (A = 1) systems can be studied numerically. Such 
studies for atomic and molecular species [12j [T4HT7] can provide useful insight into the failings of standard density- 
functional approximations and provide data for the construction and evaluation of new forms, based on the modeling 
of the adiabatic connection [TSl ESH5D] . 



Having established a convex formulation of CDFT in Sec. HIE it is possible to calculate the adiabatic connection 
in CDFT in a manner similar to that of standard DFT. Details of the adiabatic connection for CDFT have been 
presented previously by Liu [2S]. The ground-state energy functionals at interaction strength A are given by 

E x [v,A]=mfTr(TH x [v,A}), (49) 

E\[u, A] = E\[u — iA 2 , A], (50) 

where we have introduced the Hamiltonian 

k k 

From the exact functional E x [u,A] or an accurate approximation to it, the adiabatic connection may be studied in 
terms the corresponding Vignale-Rasolt functional 

^vr,a[pJ p ] = sup [E x [u, A] - (p\u) - (j p |A)] , (52) 

ti.A 

which needs to be evaluated for a fixed density (p, j p ) and different values of A in the interval < A < 1. Typically, the 



density is the ground-state density for some external potential (v cx ±, A cxt ) at A = 1. The optimization in Eq. (52) is 
trivial for A = 1 since the optimal potential is then (v cxt + ^A^ xt , A cxt ); for < A < 1, the optimization is non-trivial. 
In particular, for A = 0, the optimal potential is the Kohn-Sham potential (u s , A s ) given by 

Ms = «ext + vj + v xc + \Al = v s + \Al, (53) 
A s = A 0Xt + A xc , (54) 

where the classical Coulomb or Hartree potential wj and the exchange-correlation potential w xc are the functional 
derivatives of the corresponding energy components as in standard Kohn-Sham DFT. The exchange-correlation 
contribution to the vector potential is defined as 

A xc = (55) 

where differentiability in relevant directions is assumed. These scalar and vector potentials then enter the CDFT 
Kohn-Sham equations [26], which may be re- written in terms of (u s , A s ) as 

[\p 2 + i{p, A s } + u s ] <p p = e p ifp. (56) 



12 



If the spin dependent B(rfc) • S term is included in the Hamiltonian of Eq. (10) with the modified interactions of 



Eq. ( 48 1 , then similar arguments apply. Two-component spinors rather than one-particle orbitals then occur in the 
Kohn-Sham equations, allowing for a treatment of non-collinear magnetism. 

Even in the absence of external magnetic fields, violations of non-interacting w-representability — that is, the exis- 
tence of a ground-state density of the fully interacting Hamiltonian Hi[v,A] that cannot be reproduced by Slater- 
determinantal ground-states of the non- interacting Hamiltonian Hq[v s , A s ] — have been shown to be common in two- 
electron systems [51] ■ In general, an extended Kohn-Sham formalism, allowing for an ensemble description and 
fractional occupation numbers, is therefore required in CDFT as well as in standard DFT. 

To facilitate the optimization of i 7 vR,A[p, j P ] at a general interaction strength A, we restrict our attention to classes 
of potentials that can be parameterized in a simple way — for example, as linear combinations of basis functions. The 
most direct way to benefit from the concavity of E[u, A] is to parameterize u rather than v in the affine form 

u(r)=tw(r)+X>/*( r )» (57) 
t 

A(r)=A ref (r)+^c t5t (r). (58) 
t 

The use of u rather than v eliminates the A 2 term and the associated quadratic dependence on c t , thereby simplifying 



the equations obtained upon substitution in Eq. ( 52 1 . Importantly, it also ensures that all stationary points are 
true global maxima. To perform optimizations similar to those in Refs. [T3TU7| . all that remains is to be able to 
calculate the ground-state energy E\[u, A] with sufficient accuracy and to choose appropriate basis functions {ft} 
and {<?*}. Derivatives with respect to the expansion coefficients bt and c t may then be used in a quasi-Newton 
procedure analogous to that in Ref. 13 . The choice of basis functions and reference potentials raises the issue of a 
balanced descriptions of u and A, as well as the asymptotic limits and reduction of gauge freedom inherent in the 
finite basis set. Although important for practical implementation, these issues are beyond the scope of this article. 



I. A note on spin densities 



The formulation of CDFT with (u, A) as the basic potential and with the intrinsic and extrinsic energies expressed 
as mutual Legendre-Fenchel transforms makes the introduction of spin straightforward. The addition of the spin- 
Zeeman operator to the Hamiltonian introduces an energy term containing the spin density m(r) paired with the 
magnetic field, (m|B) = (m|V x A). A partial integration transfers the curl operator to the spin density and gives a 
surface term if the integration domain is finite. Neglect of the surface term leads to the energy term 

(m|V x A) = (V x m|A) (59) 

and a theory with (p, j m ) as the basic density. Here, 

jm = j p + V x m, (60) 

is the sum of the paramagnetic current and the spin current V x m. All results above remain valid with j p replaced 
by j m and with suitable modifications of the Hamiltonian and definitions of Vjv and B^. 

The use of j m as a basic density has been discussed by Capelle and Gross [55] as a way to translate functional 
between spin-density functional theory (SDFT) and CDFT. The alternative formulation of current-spin DFT (CSDFT) 
in terms of the charge density p, the spin density m, and the paramagnetic current density j p as separate basic 
variables is a less attractive formal theory since two independent potentials such as (u, A) cannot be conjugate to 
three independent densities (p, m, j p ). On the other hand, construction of practical approximate exchange-correlation 
functionals may be substantially more difficult with the basic variables p and j m . 



IV. THE PHYSICAL CURRENT AS A BASIC VARIABLE 



Mathematically, it is not surprising that a theory formulated in terms of magnetic vector potentials and wave 
functions — both gauge-dependent objects — makes use of gauge-dependent basic variables. Indeed, with gauge- 
dependent notions being so deeply entrenched in the theory, it is not trivial to construct a useful reformulation that 
features only gauge-invariant basic variables. On the other hand, from a physical point of view, it is somewhat 
unappealing that the paramagnetic current, rather than the physical current, arises as a basic density. Therefore, 
some authors have attempted the construction of an alternative CDFT, with the physical current as a basic variable. 
In particular, Pan and Sahni have gone far in arguing that the paramagnetic current density, in some sense, cannot 
correctly be regarded as a basic CDFT variable [22 25 . 
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A. Do physical densities determine potentials up to a gauge? 



An important question that arises in CDFT is to what extent an HK theorem is possible when the physical densities 
(p,j) are chosen as the basic densities — that is, whether (vi, Ai) ^ (uj, A 2 ) always implies that (piji) ^ (p2,h)- 

Two simple observations lend plausibility to this claim. First, for one-electron systems, (v, A) can be found explicitly 
given (p, j). Writing the one-electron wave function as 

V>(r) = i?(r)e iS(r) , (61) 

where R and S are real valued, we obtain the density p = R 2 and physical current density j = />(VS + A). Because 
of the identity V x VS = 0, the physical densities p and j determine the external magnetic field, 

B(r) = VxA(r) = Vx^|. (62) 

P( r ) 

The scalar potential v may be determined up to a constant from the eigenvalue equation ( ^tt 2 + v — E)ip = and the 
observation that 7r = — iV + A = — iV + j/p — VS, yielding 



i(-iv + M) 2 + „ W - £ 



P 1/2 (r)=0, (63) 



from which v(r) — E is uniquely determined. Hence, for a one-electron system, B and v are both determined by p 
and j. For an ./V-electron system, we obtain more generally 

Vx^=i/(r)+B(r), (64) 

where v = V x p~ 1 j p is the vorticity introduced by Vignale and Rasolt [26 . Vorticity is a gauge-invariant quantity, 
but it is not clear whether it can be uniquely reconstructed from (p,j), without a priori knowledge of the external 
magnetic field B. 

Following Pan and Sahni [55], the second observation is that, while there is no HK theorem for CDFT in terms of 
paramagnetic densities, the counterexamples to the existence of an HK theorem for paramagnetic densities (such as 



the harmonic-oscillator system in Sec. Ill C I do not preclude an HK theorem for physical densities. 

To see this, note that these counterexamples all exhibit different potentials (wi,Ai) ^ (i^Aa) with the same 
ground state tp\ = ip2 and therefore the same paramagnetic density (pi,j P i) = (p2, jp2)- However, according to the 
original HK theorem, this situation is impossible when Ai = A 2 since this would imply v 2 ^ v± + c and therefore 
Pi 7^ P%- We must therefore assume that K\ ^ A 2 in the counterexamples. However, from the assumption j pl = j p2 , 
it then follows that ji — J2 = Pi(Ai — A 2 ) ^ 0. In short, if (ux,Ai) * {v%, A 2 ) share the same ground-state wave 
function, then the physical densities are not the same. It remains to explore whether it is possible to have different 
ground states ip\ ^ ip2 but the same physical densities (pi,ji) = (p2jj2); if possible, then no HK theorem exists for 
the physical current densities. 

General arguments for an HK theorem for physical current densities have been put forth by Pan and Sahni [55J [53] 
and by Diener |44J. However, as discussed below, neither of these arguments amounts to a rigorous proof. To our 
knowledge, the existence of an HK theorem for physical current densities therefore remains open. 



B. Pan and Sahni's argument 

The standard HK theorem of DFT states that !)/u'|c implies p' ^ p' . The proof has two parts: first, it is shown 
that scalar potentials that differ by more than a constant must have different wave functions ip and ij) 1 ; second, this 
result is combined with the Rayleigh-Ritz variation principle to show that the assumptions p = p' and d^d' + c lead 
to a contradiction. It follows that » / u' + c implies p ^ p' . 

For CDFT with physical densities, a proof along the same lines has been attempted by Pan and Sahni [55]. The first 
part of their argument establishes that different potentials (v,A) (v',A') cannot yield the same physical density 
G°0)jo) = (Po)Jo) if t ne ground-state wave functions are the same. That is, without loss of generality, it may be 
assumed in an HK-type argument that the wave functions are different. The second part of their argument seeks 
to establish, by the Rayleigh-Ritz variation principle, that gauge-inequivalent potentials cannot have both different 
ground-state wave functions and the same physical ground-state density. To this end, two potentials (u, A), (y' , A') € 
Vjv, (v, A) ^ (i/,A'), are considered and it is argued that a contradiction arises from the following assumptions: 
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(a) The potentials have the same physical density (po,jo). 

(b) The potentials differ by more than a gauge transformation, (V, A') (v, A). 

No contradiction can result from (a) alone, so the assumption (b) is crucial for any reductio ad absurdum argument 
to succeed. For example, because of gauge freedom, (a) and the negation of (b) are certainly not contradictory since 
this corresponds to the perfectly consistent situation where two potentials that differ by a gauge transformation give 
rise to the same physical density. A correct proof must therefore contain at least one step that makes use of (b). 
However, while stated as an assumption, (b) is in fact never used in Pan and Sahni's argument. The argument must 
consequently be invalid. 

In more detail, we here identify an erroneous step in the reasoning in Ref. in particular Eqs. (35)-(40). With 
ijj and ip' denoting the ground states of the two potentials, we obtain Eq. (35) from Ref. [52]: 

E = {i)\H[v,A]\4>) < (ip'\H[v,A}\ip'). (65) 

The identity H[v, A] = H[v', A'] + (H[v, A] — H[v', A']) yields for the expectation value on the right-hand side 

W\H[v, AM) =E' + (V/|i7r 2 + v - \-k' 2 - v'W) 

= E' + ($'\\{-n,A}+v- \A 2 - \{k',A'}-v' + \A' 2 W). 

It is here important to distinguish between tv = p + A and 7r' = p + A' since the representation of the mechanical 
momentum operator is gauge and vector-potential dependent. To proceed, we explicitly write out the corresponding 
physical current density operators, which for the two potentials are given by 



j(r) = \ X> fc *(r fc - r) + 6(r k - r)n k ), (67) 

k 

j'(r) = \ J2«S(r k - r) + S(r k - r)n' k ). (68) 

k 

From these, we may calculate the physical current (assumed to be the same in the two cases) and its interaction with 
some vector potential a(r) as 

j„(r) = (^|j(r)|V) = (^li'W) (69) 
(j |a) = Mlfr.aM = (70) 

where it is important to use primed or unprimed quantities consistently. To ensure that we are using the correct 



current-density operator for ip' in Eq. (66), we insert the identity {tv, A} = {tv 1 + A — A', A} yielding, 



E < E' + ( Po \Av) + (j 1 A A) + l( Po \AA 2 ). (71) 



where we have introduced Av = v — v' and A A = A — A' and also used the inequality in Eq. ( 65 ) 



Carrying out the above argument with primed and unprimed variables interchanged, we obtain the strict inequality 

E + E' <E' + E +{p Q \AA 2 ), (72) 

where the last term does not vanish since AA 2 is symmetric in the primed and unprimed variables. In agreement 
with the HK theorem, a contradiction arises if (and only if) A = A'. When A(r) and A'(r) differ where p(r) ^ (as 
a result of different physical fields or by a gauge transformation), no contradiction arises. 

In the argument given in Ref. |22j . the authors incorrectly identify (ip' \j(r)\ip') and (tp' \j' (r)\ip') , leading to their 



Eq. (38), which differs from Eq. (71) above by the replacement of ^(po\AA 2 ) with ^(po\A' 2 — A 2 ). Since their term 
is antisymmetric rather than symmetric in the primed and unprimed variables, the authors obtain a contradiction 
E + E' < E + E' irrespective of A and A', leading to the unjustified conclusion that (poi jo) determines (v, A). 



C. Diener's argument 



An earlier attempt to prove an HK-type theorem for physical currents by Diener [33] invokes an intriguing strategy 
for eliminating the term (pq\AA 2 ) from Eq. (72). The key idea is to replace the external vector potential A by an 
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effective vector potential a c ff[j,-0] defined to reproduce a given physical current density j from a given wave function 
if] so that 

j = j P ;V + /=tya e ff [j, # (73) 
Consider now the following rearrangement of the kinetic energy term 

tt 2 =p 2 + {p,A} + A 2 

= p 2 - a 2 cS + {p + a cff , A} + (a cff - A) 2 . (74) 

By definition, ^ (-0|{p + a e g[j, A}|-0) = (j|A). For a prescribed physical current density j , the expectation value 
of an arbitrary wave function is then 

(i,\H[v, AM = (^(p 2 - a 2 cS ) + W\j>) 

+ (jo|A) + (p| i ;) + i(p|(a cff -A) 2 ). (75) 

The total current evaluated using the effective momentum operator p + a c ff instead of the true mechanical momentum 
operator p + A always evaluates to the prescribed current jo- Diener then defines a universal density functional of 
the form 



Pb[pJ] = mf (mMi^M, 

By exploiting this functional, Diener derives the inequality 

(^ |#[«,A]|V>o) < W \H[v,A]\il/ ) - \(po\AA 2 ), 



(76) 
(77) 

(78) 



in lieu of the usual strict Rayleigh-Ritz inequality underlying the HK proof as given in Eq. ( 65 ) . When is the 
ground state of potentials that differ by more than a gauge, (i/,A') oo (w,A), the above non-strict inequality is 
(without further ado) taken to be strict in Diener's presentation. If a strict inequality is accepted, a standard reductio 
ad absurdum proof is possible because the last term of Eq. ( |72| ) cancels to yield the contradiction E + E' < E' + E. 

We now consider two technical problems not addressed by Diener. The first is that the expectation value being 
minimized in Fd is not bounded from below for densities that vanish at some point. Although unusual, such ground- 
state densities can arise for small molecules in strong magnetic fields. Consider a wave function <f>o giving rise to 
a density po that vanishes at some point in space O. Let us also introduce spherical coordinates (r, 9, <f>) about O. 
Then, for the wave functions 



$ m (x 1; . . . , XJV ) = e Mm ^ 0fc $o(xi, . . . ,x w ) 



(79) 



with integer m, we see that $ m and $0 gi ye r i se to the same density po but to different paramagnetic current densities 
related by 



Jp;m 



PO + jp;0, 



r sm(8) 

where <f> is the unit vector in the direction specified by <fi. As a result, 

m4> 



a m = a cff [j, $ m ] = a cff [j, $ 



r sm(6) 



a 



m<j) 
■ sin(0) ' 



(80) 



(81) 



Calculating the expectation value of the effective Hamiltonian in Eq. (77), some terms arising from cancel terms 
arising from |a^, leaving 



($o\Hp 2 



f j 

r sm(0) 



dr. 



(82) 
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For physical currents j with a nonzero last integral, for example a circular current, the expectation value can be 
decreased without bound, demonstrating that -Fb[p,j] is not well defined on the full domain Bn- 



The second technical problem in the derivation of the inequality in Eq. (78) concerns the strategy of choosing 
a prescribed jo such that the effective vector potential equals the external vector potential. This amounts to a 
self-consistency condition [Eq. (11) of Ref . |4"4"] . 

jo = jp;* D + Poa c ff [jo, 'I'd] = j p; *[u,A] + A)A, (83) 
where <E>d minimizes ($D|^cff[jo, 3>d]|$d) and ^[t^Alis the ground state of the external potentials giving rise to 



(po, jo). Diener's derivation of the inequality in Eq. (78) hinges on the ability to reproduce the paramagnetic current 
density from a ground state \&[u,A] of H[v,A] by a minimizing wave function in i<b[pojjo]- However, it is unclear 
whether this is always possible. 

D. Constrained search with physical currents 

In lieu of a rigorous proof for an HK-type theorem for physical currents, it is possible to proceed by conjecturing 
such a result and exploring its consequences. Additionally, an important point is that a mapping from ground state 
densities to potentials may not be required for a formulation of CDFT, provided that the theory can be constructed 
by some other means such as a constrained search or Legendre-Fenchel transformation formalism. We shall here 
explore such issues in the framework introduced by Pan and Sahni [22 25 in more detail. 

A complication due to the choice of the physical current as a basic variable is that constraints of the type "ip — > p, j" 
require explicit reference to a vector potential, because the physical current is not determined by the wave function 
alone. Care must therefore be exercised when developing a constrained-search formalism for physical currents. 

For each magnetic field under consideration, we fix a gauge. Hence, we choose a mapping 

B(r)h->a[B](r) (84) 

from magnetic fields to magnetic vector potentials, and also fix the constant shift of scalar potentia ls in so me way. 
From the conjecture that a ground state density (/?, j) uniquely determines a gauge class of potentials (v, A), we may 
now write mappings 

(p,j)o(V«,B)o(«,a[B]) (85) 

Hence, we may write the scalar potential, the external magnetic field and its vector potential as functionals v = w[p, j], 
B = b[/o,j], and A = a[b[p,j]] of the physical densities; these are representatives of the equivalence class (v, A). 

Given potentials (v, A), we may determine the corresponding ground state ipo an d physical ground-state densities 
(p,j) € Bn, and express the energy as 

E[v, A] = (ij \y 2 + W\ifo) + (p\v - \A 2 ) + (j|A). (86) 

Using the mapping from densities to potentials, a universal density functional could in principle internally reconstruct 
the ground state potentials, and obtain the intrinsic energy via 

F PS [p,j] =£[Mp,j],a[b[p,j]]] - (p\w[p,j] - §a[b[pj]] 2 ) 

-(j|a[b[p,j]]). 1 ) 

Furthermore, we may define a non-universal constrained-search functional that depends explicitly on the external 
vector potential, 

F£f[p,3,A}= inf Tr(r(ip 2 + W)) 

rn-(/3j-pA) /gg\ 

= Fy R [p,j- pA], 

where we note that j — pA — j p . In principle, given the conjectured one-to-one mapping between potentials and 
physical densities, a functional of physical densities may internally reconstruct the corresponding potentials. This 
observation is relied on in the definition of Fps and it is tempting to consider also a universal constrained-search 
functional that exploits this idea, 

i^ s 2 [p,j]= .inf Tr(r(ip 2 + M0) 

!>>(pj-pa[b[p,j]]) ^gg-j 

= F VR [p,j - pafb^J]]]. 
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Consider now the reformulation of the minimization in Eq. (Ill as a nested minimization over physical densities 
and wave functions. 



F£lVJ',A] 



E[v,A] = mi inf Tr(r(|p 2 + W)) 



1 

2 J 



+ ( P >-iA 2 ) + (j'|A) 



(90) 



A constrained-search-like expression is thus possible, but the functional that arises from the nested minimization is 
the non-universal functional -Fpg 1 . A minimization over the universal functional -Fp| 2 is here not equivalent, 



V,j' 



E[v,A] ^ inf 



r^p'J'-p'a[b[p',j']] 

+ {p'\v-\A 2 ) + {i'\A) 



inf Tr(r(ip 2 + W)) 



(91) 



To establish the non-equality, consider the introduction of a new scalar potential u_ = v — \A 2 . Written in terms of 
it_, the right-hand side of Eq. (91 1 has the form of a Legendre-Fenchel transformation and consequently a functional 
of (w_, A) that is concave by construction. However, writing the harmonic oscillator energy of Sec. HIE in terms of 



this scalar potential, we see that the resulting energy functional is not concave in (u_ , A) (it is in fact only concave 
in scalar potentials u\ = v + XA 2 for A > 5). Hence this non-concave energy functional cannot equal the concave 
functional defined in terms of -Fp| 2 . In previous work (see, e.g., Eqs. (24)-(28) of Ref. [IS]), these two non- equivalent 
functionals were conflated. 

The constrained-search approach to CDFT with the physical current as a basic variable is substantially complicated 
by the fact that a wave function does not determine the physical current. It does not appear possible to construct 
a functional that is both universal and admits a straightforward variation principle. Moreover, unlike the situation 
when j p is a basic variable, a simple redefinition of the scalar potential does not simultaneously yield a concave energy 
functional and linear pairing between potentials and physical densities. An avenue for further study could be to 
consider the properties of a universal intrinsic energy functional of the physical densities. Such a functional cannot 
be convex in the pair (p, j) as a whole, but could possibly be convex in p (for a fixed j) and concave in j (for a fixed 
p), which would enable at least partial Legendre-Fenchel transformations. 



V. CONCLUSIONS 



In this work we have compared formulations of CDFT based on different choices of basic variables. While the usual 
formulation in terms of (p,j P ) and {v,A) does not lead to a one-to-one mapping between potentials and densities, 
it allows for the construction of a useful density functional theory via a constrained search approach. A drawback 
with this choice of basic variables is that E\v, A] is not concave, making a full Legendre-Fenchel transform treatment 
infeasible. However, as shown here, concavity can be restored by considering the conjugate variables (p, j p ) and (u, A), 
where u = v + \A 2 . This allows the application of convex analysis in analogy with Lieb's formulation of standard 
DFT [3 . Such a formulation is particularly natural in the context of the study of adiabatic connections for CDFT 
functional construction, which can be done in a manner similar to that undertaken for standard DFT [321 IT4HT7] . 
The information garnered from such analysis would be suitable for comparison with Kohn Sham implementations of 
CDFT based on functionals of (p,j P ). 

Alternative CDFT formulations based on the physical densities (p,j) have also been critically examined. Pan and 
Sahni's 22-25] recent attempt to formulate such a theory is found to be unsatisfactory. Both their attempt to prove 
the existence of a one-to-one mapping between potentials and densities and their constrained-search formulation were 
found to be flawed. An earlier attempt to formulate a CDFT in terms of the physical current by Diener [23] has also 
been examined, and technical difficulties with this approach have been highlighted. Despite the appealing physical 
motivation behind this choice of basic variables, we thus find that a formal justification for such a framework is 
currently lacking. Furthermore, while it remains open whether or not an analogue of the Hohenberg-Kohn theorem 
holds for the physical current, other aspects of standard DFT such as the variation principle, the constrained-search 
formalism, and formulations in terms of Legendre-Fenchel transformations do not straightforwardly carry over to this 
type of CDFT. We conclude that the most common formulation in terms of (p, j p ) is presently the most convenient 
and viable formulation of CDFT. 
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APPENDIX 



Let X be a normed vector space and X* its dual — that is, the the linear space of all continuous linear functionals 
on X. A function / : X — > E U {+00} is said to be convex if it satisfies the relation 

f(\ Xl + (1 - X)x{) < A/On) + (1 - X)f(x 2 ) (92) 

for all < A < 1. The effective domain, dom(/), is the set of x 6 X for which f(x) < +00. A function / is lower 
semi-continuous at xq if, for any e > 0, there exists S > such that f(x) > f(xo) — e whenever \x — xq\ < S; / is lower 
semi-continuous if it is lower semi-continuous at all x € X. 

A lower semi-continuous convex function / : X — >• M U {+00} may be represented by its conjugate /* : X* —> 
M. U {+00} in the manner 

f*(y) = sup[(x\y) -/(*)], (93) 

/(s)= supply) (94) 

sex* 

where /* is also lower semi-continuous and convex. The dual function yo g A* is said to be a subgradient of / at 
x G X if it satisfies the inequality 

f(x)>f{x ) + (x-x \y ), Viel. (95) 

The set of all subgradients of / at xq is called the subdifferential df(xo) and is a (possibly empty) convex subset of 
A*. Subgradients and subdifferentials of /* are defined in an analogous manner. The function / and its conjugate 
satisfy Fenchel's inequality f(x) + f*(y) > {x\y), which is sharpened into the equality f(x) + f*(y) = (x\y) whenever 
the equivalent reciprocal relations 

yedf(x) ^ xedf*(y) (96) 

are satisfied. 

A function / :I->1U {— 00} is said to be concave if — / is convex. Its conjugate is defined in the same manner 
as for convex functions but with sup replaced by inf; likewise, subgradients and subdifferentials are defined as for a 
convex function but with the inequality sign reversed in Eq. ( 95 1 . 
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